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Abstract. Consider a planar, bounded, m-connected domain f2 and let dil be its boundary. 
Let T denote a cellular decomposition of 17 U dfl where each 2-cell is either a triangle or a 
quadrilateral. We construct a new decomposition of Q U d£l into 1Z, a finite, disjoint union 
of quadrilaterals. The construction is based on utilizing a pair of functions on and 
properties of their level curves. The first function is obtained as the solution of a Dirichlet 
boundary value problem defined on . The second function is obtained by an integration 
scheme along the level curves of the first, and is called the conjugate function. It turns out 
that the pair of functions can be made orthogonal and harmonic in a combinatorial sense. 
As an application, we advance beyond the main results in |18j . 



0. Introduction 

0.1. Perspective. In his attempts to prove uniformization Riemann suggested considering 
a planar annulus as made of a uniform conducting metal plate. When one applies voltage to 
the plate, keeping one boundary component at voltage k and the other at voltage 0, electrical 
current will flow through the annulus. The equipotential lines form a family of disjoint simple 
closed curves foliating the annulus and separating the boundary curves. The current flow 
lines consist of simple disjoint arcs connecting the boundary components, and they foliate the 
annulus as well. Together the two families provide "rectangular" coordinates on the annulus 
that turn it into a right circular cylinder or a (conformally equivalent) circular concentric 
annulus. 

In [H], Brooks, Smith, Stone and Tutte studied square tilings of rectangles. They defined 
a correspondence between square tilings of rectangles and planar multigraphs endowed with 
two poles, a source and a sink. They viewed the multigraph as a network of resistors in which 
current is flowing. In their construction, a vertex corresponds to a connected component of 
the union of the horizontal edges of the squares in the tiling; an edge appears between 
two such vertices for each square whose horizontal edges lie in the corresponding connected 
components. Their approach is based on Kirchhoff's circuit laws which are widely used 
in the field of electrical engineering. It is interesting to recall that it was Dehn [Hj, who 
was the first to show (in 1903) a relation between square tiling and electrical networks. 
In [22], combining a mixture of geometry and probability, Kenyon used the fact that a 
resistor network corresponds to a reversible Markov chain to show a correspondence between 
a class of planar non- reversible Markov chains and trapezoid tilings. In [TS] and [19] , we 
addressed (using methods that transcend Dehn's ideas) the case where the domain has higher 
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connectivity. This provides a first step towards an approximation of conformal maps from 
such domains onto a certain class of surfaces. 

In [27] , starting with a triangulation of a rectangle, Schramm employed discrete extremal 
length arguments and constructed a tiling of a rectangle by squares that preserves the contact 
graph of the given triangulation. In [TTJ, Cannon, Floyd and Parry used extremal length 
arguments (similar to those in [27]) and gave a proof of a similar theorem. Another proof 
of a slight generalization of [27] was given by Benjamini and Schramm [7] (see also [8] for 
a related study). Both [27] and [11] are widely known as the "Finite Riemann Mapping 
Theorem," which serves as the first step in Cannon's "Combinatorial Riemann Mapping 
Theorem" (see |10j). We will comment about the relation of the above mentioned works to 
Theorem 10.41 in 13.1.11 

One may view the above mentioned theorems by Schramm, and by Cannon, Floyd, and 
Parry as attempts to provide various types of discrete versions of the conformal map obtained 
by the Riemann Mapping Theorem. 

0.2. Combinatorial rectangular coordinates. In this paper, we will follow Riemann's 
perspective on uniformization by constructing "rectangular" coordinates from given com- 
binatorial data. The foundational modern theory of boundary value problems on graphs 
enables us to provide a unified framework to the theorems mentioned above, as well as to 
more general situations. The important work of Bendito, Carmona and Encinas (see for 
instance [I], [5] and [6]) is essential for our applications, and parts of it were utilized quite 
frequently in [17], [18], [19], this paper, and its sequel [20] . 

Consider a planar, bounded, m-connected domain Q and let dQ be its boundary. Hence- 
forth, let T denote a cellular decomposition of Q U dQ where each 2-cell is either a triangle 
or a quadrilateral. We will construct a new decomposition of QUdQ into 71, a finite, disjoint 
union of quadrilaterals. We will show that the space of quadrilaterals can be endowed with a 
finite measure thought of as a combinatorial analogue of the Euclidean planar area measure. 

Next, we construct a pair (S^, /) where Sq is a special type of a genus singular flat 
surface having m boundary components, which is tiled by rectangles and is endowed with 
/i, the canonical area measure induced by the singular flat structure. The map / is a 
homeomorphism from (Q, dQ) onto Sq. Furthermore, each quadrilateral is mapped to a 
single rectangle and its measure is preserved. 

The proof that / is a homeomorphism as well as the construction of a measure on the 
space of quadrilaterals depend in a crucial way on constructing a pair of harmonic functions 
on TZ^ and a few properties of their level curves. 

Let g denote the solution of a Dirichlet boundary value problem defined on T^°' (see 
Definition II .4p . We will start by extending g to the interior of the domain: affinely over 
edges in and over triangles T^ 2 \ and bilinearly over quadrilaterals in T^. We will 
often abuse notation and will not distinguish between a function defined on and its 
extension over \T\. 

For the applications of this paper, particularly in creating "rectangular" coordinates in a 
combinatorial and topological sense, it is necessary to define a new function 9 on . This 
function will be defined on an annulus minus a slit (or a quadrilateral) and will be called 
the conjugate function. It is obtained by integrating the normal derivative of g along its 
level curves (Definition II. ip . The normal derivative of g is initially defined only at vertices 
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that belong to dQ. The simple topological structure of the level curves of g permits the 
extension of the normal derivative to the interior, and thereafter its integration. In the case 
of an annulus, these level curves are simple, piecewise linear closed curves that separate the 
boundary curves and foliate the annulus. Definition 12.61 will formalize this discussion. 

The analysis of the level curves of 9 is the subject of Proposition 12.211 Their interaction 
with the level curves of g is described in Proposition 12. 221 In the case of an annulus, the level 
curves of g form a piecewise-linear analogue of the level curves of the function u(r, <fi) = r, 
and those of 9 form a piecewise-linear analogue of the level curves of the function v(r, <fi) = <fi. 

We now turn to a more technical description of this paper. 

A simple quadrilateral is a closed topological disk with precisely four vertices on its bound- 
ary. For h : — > ffi. recall that dh(e) = h(e + ) — h(e~) for any oriented edge e = [e~,e + ]. 
A slit in an annulus is a fixed, simple, combinatorial path in along which g is monotone 
increasing which joins the two boundary components (Definition 12.11) . 

We now make 

Definition 0.1 (A pair of combinatorial orthogonal functions). Let (Q,dQ,T) be given, 
where Q is an annulus minus a slit. A pair of non- negative functions <fi and ip defined on 
7~(°) will be called combinatorially orthogonal, if there exist a cellular decomposition 1Z of 
U <9f2, where each 2-cell is a simple quadrilateral, and extensions of <p and ij) to 1Z^\ such 



Such a decomposition will be called a rectangular combinatorial net. 

We record the properties of the pair which are essential for the applications of this paper. 

Theorem 0.3. The pair {g, 9} is combinatorially orthogonal. Furthermore, there is a choice 
of conductance constants on such that both functions are combinatorially harmonic. 

0.3. Geometric models and maps. We now turn to describing some applications of the 
pair {g, 9} and the rectangular net TZ. In the course of the proofs of our main theorems, we 
will first construct a new decomposition of Q into a rectangular net, TZ, then a model surface 
which is, when m > 2, a singular flat surface tiled by rectangles. Finally we will construct a 
map between the domain and the model surface and describe its properties.. 

Let us start with the fundamental case, an annulus. Given two positive real numbers 
7*1 and ?"2 and two angles <f)i,4>2 £ [0, 27r), a bounded domain in the complex plane whose 
boundary is determined by the two circles u(r, <$>) = r\ and u(r, <fi) = 7*2 and the two radial 
curves v(r,<f)) = <fti and v(r,<f>) = 4>2, will be called an annular shell. Let \i denote Lebesgue 
measure in the plane. In the statement of the next theorem, the measure v which is described 
in Definition 13.11 is determined by g and 9. The quantity period^) is an invariant of 
9 which encapsulates integration of the normal derivative of g along its level curves (see 
Definition I2HSD- 

Our first main theorem is 



that 



(0.2) 




Theorem 0.4 (A Dirichlet problem on an annulus). Let A be a planar annulus endowed 
with a cellular decomposition T , and let dA — E\ U E^. Let k be a positive constant and let 
g be the solution of the Dirichlet boundary value problem defined on (A,dA,T). 
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Let S A be the concentric Euclidean annulus with its inner and outer radii satisfying 
(0.5) {n,r 2 } = {l,27rLength(L(^))} = {l,27rexp ( — j^— k)}. 

Then there exist 

(1) a tiling T of S A by annular shells, 

(2) a homeomorphism 

f:(A,dA,K)^(S A ,dS A ,T), 

such that f maps each quadrilateral in VS 2 ^ onto a single annular shell in S A , f 
preserves the measure of each quadrilateral, i.e., 

u(R) = n(f(R)), for all R G K {2 \ 

and f is boundary preserving. 

The dimensions of each annular shell in the tiling are determined by the boundary value 
problem (in a way that will be described later). In our setting, boundary preserving means 
that the annular shell associated to a quadrilateral in 1Z with an edge on dVL will have an 
edge on a corresponding boundary component of S A . 

Our second main theorem is Theorem 14.41 It provides a geometric model for the case 
m > 2. The model surface that generalizes the concentric annulus in the previous theorem 
first appeared in [18]. It is a singular flat, genus zero compact surface with m > 2 boundary 
components with conical singularities. It is called a ladder of singular pairs of pants. 

We obtain a topological decomposition of Q into simpler objects, annuli and generalized 
annuli, for which the previous theorem and a slight generalization of it may be applied. The 
second step of the proof is geometric. We show that it is possible to glue the pieces along 
common boundaries in a length preserving way. This entails a new notion of length which 
will be introduced in Definition 12.361 

0.4. Advances beyond [IS] . From [18] we use the description of the topological properties 
of singular level curves of the Dirichlet boundary value problem. The most significant one is 
a description of the topological structure of the connected components of any singular level 
curve of the solution. A study of the topology and geometry of the associated level curves 
and their complements is carried out in [IS1 Section 2]. 

The advances of this paper beyond [18] are due to the fact that we are able to make a 
change of charts on Q. That is, we view the given cellular decomposition T as a set of initial 
charts. Any component of the complement of a singular level curve is either an annulus or 
an annulus with one singular boundary component. We obtain a rectangular net on it using 
level curves of {g, 9}. 

Once this step is accomplished, we are able to advance beyond the main results of [18] as 
well as the other works mentioned in Subsection lO.il We will show that the map we construct 
from the domain to a model surface is a homeomorphism onto its image, and furthermore 
that it preserves a combinatorial two-dimensional measure. 

Recall that in the theorems proved in [IS] (as well as in [12]), the analogous mapping to / 
was proved to be an energy-preserving map (in a discrete sense) from onto a particular 
singular flat surface. It is not possible to extend that map to a homeomorphism defined on 
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the domain. Furthermore, the natural invariant measure considered there is one- dimensional 
(being concentrated on edges). 

In [20], we will continue our work by showing that a scheme of refining the cellular de- 
composition, coupled with a particular choice of a conductance function in each step, leads 
to convergence of the mappings constructed in each step, to a canonical conformal mapping 
from the domain onto a model surface. 

0.5. Organization of the paper. A modest familiarity with [18] will be useful for reading 
this paper. For the purpose of making this paper self-contained, a few basic definitions and 
some notations are recalled in Section HJ and results from [TH] are quoted as needed. In 
Section [2J the main technical tools of this paper, a conjugate function and a rectangular net, 
are constructed on an annulus minus a slit. In Section [3j the cases of an annulus and an 
annulus with one singular boundary component are treated, respectively, by Theorem 10.41 
and Proposition 13.221 Due to the reasons we mentioned above, these are foundational for the 
applications of this paper and of [20] as well. Section[4]is devoted to the proof of Theorem 14.41 

Convention. In this paper we will assume that a fixed affine structure is imposed on 
(Q,dQ,T). The existence of such a structure is obtained by using normal coordinates on 
(Q, dfl, T) (see [2H1 Theorem 5-7]). Since our methods depend on the combinatorics of the 
triangulation, the actual chosen affine structure is not important. 

Acknowledgement. It is a pleasure to thank Ted Shifrin and Robert Varley for enjoyable 
and inspiring discussions related to the subject of this paper. 



In this section we briefly review classical notions from harmonic analysis on graphs through 
the framework of finite networks. We then describe a procedure to modify a given boundary 
problem and T. The reader who is familiar with [IS] or [TH] may skip to the next section. 

1.1. Finite networks. In this paragraph we will mostly be using the notation of Section 
2 in [3]. Let T = (V,E,c) be a planar finite network, that is, a planar, simple, and finite 
connected graph with vertex set V and edge set E, where each edge (x, y) G E is assigned a 
conductance c(x,y) = c(y,x) > 0. Let V(V) denote the set of non-negative functions on V. 
Given F C V we denote by F c its complement in V. Set V(F) = {u G V(V) : S(u) C F}, 
where S(u) = {x G V : u(x) ^ 0}. The set 5F = {x G F c : (x,y) G E for some y G F} is 
called the vertex boundary of F. Let F = F[J5F and let E = {(x, y) G E : x G F}. Let 
r(F) = (F, E, c) be the network such that c is the restriction of c to E. We write x ~ y if 



The following operators are discrete analogues of classical notions in continuous potential 
theory (see for instance [16] and [13]). 

Definition 1.1. Let u G V(F). Then for x G F, the function 



1. Finite networks and boundary value problems 



{x,y) G E. 



(1.2) 
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is called the Laplacian of u at x. For x G 5(F), let {2/1, 2/2, • • • , Um} £ F be its neighbors 
enumerated clockwise. The normal derivative of u at a point x G <5F with respect to a set 
F is 

(1.3) ^-(F)(x) = £ c(ar,j/)(«(aO-«(j/)). 

Finally, m G P(-F) is called harmonic in F C V if Aw(x) = 0, for all x E F. 

1.2. Boundary value problems on graphs. Consider a planar, bounded, m-connected 
region f2, and let dfl be its boundary (m > 1). Let T be a cellular decomposition of Q U dQ, 
where each cell is either a triangle or a quadrilateral. Let dQ = E\ U E2, where F\ is the 
outermost component of <9f2. Invoke a conductance function C on thus making it a 

finite network, and use it to define the Laplacian on . 

We will now describe the type of boundary value problem defined on that will be 
employed in this paper. 

Definition 1.4. The Dirichlet boundary value problem (D-BVP) is determined by requiring 
that 

(1) g\ El = k, g\ E2 = 0, and 

(2) Ag = at every interior vertex of . 
These data will be called Dirichlet data for Q. 




Figure 1.5. An example of a possible splitting scheme. 



A fundamental property which we often will use is the discrete maximum-minimum prin- 
ciple, asserting that if u is harmonic on V C V, where V is a connected subset of vertices 
having a connected interior, then u attains its maximum and minimum on the boundary of 
V (see [281 Theorem 1.35]). 

The following proposition (cf. [3, Prop. 3.1]) establishes a discrete version of the first 
classical Green identity. It played an important role in the proofs of the main theorems in 
[T7J [TS], and it also plays an important role in this paper and in its sequel [20] . 

Proposition 1.6 (The first Green identity). Let F C V and u, v G V(F). Then we have 
that 

(1.7) ^2 c ( x ,y)( u ( x ) - u (y))( v i x ) - v (v)) = \ Au(x)v(x)+ j —(f)(x)v(x). 

(x,y)eE 



x&F Jxg6(f) 9n 
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1.3. Piecewise-linear modifications of a boundary value problem. We will often need 
to modify a given cellular decomposition, and thereafter to modify the initial boundary value 
problem. The need to do this is twofold. First assume, for example, that L is a fixed, simple, 
closed level curve of the initial boundary value problem. Since L D is not (generically) 
a subset of Definition 12.361 may not be employed directly to provide a notion of length 
to L. Therefore, we will add vertices and edges according the the following procedure. Such 
new vertices will be called type I vertices. 

Let Oi, O2 be the two distinct connected components of the complement of L in Q, with L 
being the boundary of both (these properties follow by employing the Jordan curve theorem). 
We will call 0\ an interior domain if all the vertices which belong to it have g- values that are 
smaller than the g-value of L. The other domain will be called the exterior domain. Note 
that by the maximum principle, one of (D\,Oi must have all of its vertices with g- values 
smaller than the g- value of L. 

Let e G T^~> and assume that x = e D L is a vertex of type I. Thus, two new edges (x, v) 
and (u, x) are created. We may assume that v G 0\ and u G O2. Next, define conductance 
constants c(v,x) and c(x,u) by 



(1.8) 5(t,,s) = and ~ cM = c(v,u)(g(u)-g(v)) _ 

g(v)-g(x) g{u)-g{x) 

By adding to T all such new vertices and edges, as well as the piecewise arcs of L deter- 
mined by the new vertices, we obtain two cellular decompositions, To 1 of 0\ and 7e> 2 of Oi- 
Two conductance functions, Cq x and Cq 2 are now defined on the one-skeleton of these cellular 
decompositions, by modifying according to Equation fll.Sp the conductance constants that 
were used in the Dirichlet data for g (i.e., changes are occurring only on new edges, and on 
L the conductance is defined to be identically zero). One then defines (see [TSJ Definition 
2.7]) a natural modification of the given boundary value problem, the solution of which is 
easy to control by using the existence and uniqueness theorems in [3J. In particular, it is 
equal to the restriction of g to O i: for 2 = 1,2. 

Another technical point which motivates the modification described above will manifest 
in Subsection 12.31 Proposition 11.61 will be frequently used in this paper, and it may not 
be directly applied to a modified cellular decomposition and the modified boundary value 
problem defined on it. Formally, in order for Proposition 11.61 to be applied, the modified 
graph of the network needs to have its vertex boundary components separated enough in 
terms of the combinatorial distance. Whenever necessary, we will add new vertices along 
edges and change the conductance constants along new edges in such a way, that the solution 
of the modified boundary value problem will still be harmonic at each new vertex and will 
preserve the values of the solution of the initial boundary value problem at the two vertices 
along the original edge. Such vertices will called type II vertices. 

Formally, once such changes occur, a new Dirichlet boundary value problem is defined. 
The existence and uniqueness of the solution of a Dirichlet boundary value problem (see [3]) 
allow us to abuse notation and keep denoting the new solution by g. 
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2. CONSTRUCTING A COMBINATORIAL NET AND HARMONIC CONJUGATE FUNCTION ON 

AN ANNULUS WITH A SLIT 

This section has three subsections. The first subsection contains the construction of the 
conjugate function to the initial Dirichlet boundary value problem defined on an annulus. 
The second subsection supplies the proof of Theorem 10.31 Finally, in the third subsection, 
we will define the pair-flux metric and its induced length. These notions will be essential to 
the proof of Theorem 14.41 in which gluing two components of the complement of a singular 
level curve of the solution occurs. 

2.1. Constructing the conjugate function. In this subsection, we will construct a func- 
tion which is conjugate (in a combinatorial sense) to g, the solution of a Dirichlet boundary 
value problem defined on an annulus. The conjugate function will be single valued on the 
annulus minus a chosen slit. 

With the notation of the previous section and the introduction, let (A, dQ = EiUE 2 , T) be 
an annulus endowed with a cellular decomposition in which each 2-cell is either a triangle or a 
quadrilateral. Let k be a positive constant, and let g be the solution of a Dirichlet boundary 
value problem as described in Definition 11.41 All the level curves of g are piecewise, simple 
closed curves separating E\ and E 2 (see Lemma 2.8 in [18] for the analysis in the case of 
higher connectivity) which foliate A. 

Before providing the definition of the conjugate function to g, we need to make a choice 
of a piecewise linear path in A. 

Definition 2.1. Let slit (^4) denote a fixed, simple combinatorial path in which joins 
Ei to E 2 . Furthermore, we require that the restriction of the solution of the boundary value 
problem to it is monotone increasing. 

Remark 2.2. The existence of such a path is is guaranteed by the maximum principle. 
Let 

(2.3) C={L(v ),...,L(v k )} 

be the collection of level curves of g that contain all the vertices in arranged according 
to increasing values of g. It follows from Definition 11.41 that L(v Q ) = E 2 and L(vk) = E\. 

In order to obtain a single valued function, we will start with a preliminary case. To this 
end, let Q s i it denote the quadrilateral obtained by cutting open A along slit (.4) and having 
two copies of slit (.4) attached. Since in A orientation is well defined, we will denote one of 
the two copies by <9<2basc an d the other by dQ t0 p- In other words, from the point of view 
of A, points on slit (A) may be endowed with two labels, recording whether they are the 
starting point of a level curve (with winding number equals to one) or its endpoint. We keep 
the values of g at the vertices unchanged. Thus, corresponding vertices in <9<2base and dQ top 
have identical g-values. By abuse of notation, we will keep denoting by the 0-skeleton 

Of Qslit- 

For v G A \ E 2 , which is in or a vertex of type I, let L(y m ) denote the unique level 
curve of g which contains v. Let Q v be the piecewise-linear quadrilateral whose boundary is 
defined by slit (.4), L(v m ), L(v m -i) and the leftmost edge emanating from v with its endpoint 
on L(v m -i). 
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Figure 2.4. An example of a quadrilateral Q 



For v G D E 2 recall that E 2 = L(vq) is the unique level curve of g which contains v. 
Let Q v be the piecewise-linear quadrilateral whose boundary is defined by slit (.A), E 2 , L(v\) 
and the rightmost edge emanating from v with its endpoint on L{v\). 

Remark 2.5. Note that while Q v was defined in A, it is a well defined object in Q s n t as well. 

We are now ready to present 

Definition 2.6 (The Conjugate function). Let v be a vertex in or a vertex of type I. 
Let 

(2.7) tt(v) = L(v) n <9Q basc . 

We define 0(f), the conjugate function of g, as follows. 
First case: Suppose that v £ E 2 . Then 

(2.8) 9{v) = f ^(Q v )(u), 

Jn(v) on 

where the integration is carried along L(v) in the counter-clockwise direction. 
Second case: Suppose that v G E 2 . Then we define 9{v) by 

(2-9) 6{v)= f" \^-(Q v )(u)\, 

where the integration is carried along E 2 in the counter-clockwise direction. For a point 
z G Q S Ht which is not a vertex as above, 6(z) is defined by extending 9 affinely over edges 
and triangles, and bi-linearly over quadrilaterals. 
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Remark 2.10. The reasons for using Q(v), the "rightmost" and the absolute value of the 
normal derivative of g in the second case, are due to the orientation along E 2 and the 
maximum principle, respectively. Furthermore, the continuity of 9 from the right on E 2 
follows from similar arguments to those appearing in the proof of Proposition 12.131 

Remark 2.11. It is worth mentioning that the modification of the definition to the case of 
a quadrilateral, Q, is immediate. One obsereves that once a Dirichlet-Neumann boundary 
value problem is imposed (see [IE]), the level curves of the solution are disjoint, piecewise- 
linear simple curves that foliate the quadrilateral and join its bottom to its top (see [T9| 
Proposition 2.1]). Hence, the definition above would apply once we define slit ( Q) to be the 
bottom of <2 (<9Qbottom), an d then the analogue to Q v is defined as above. 

Remark 2.12. Note that we abuse notation and keep denoting by E\ and by E 2 the two 
boundary components of Q s u t , which naturally correspond to their counterparts in A. 

We now turn to studying the topological structure of the level curves of 9. 

By definition, <9Qbasc is the level curve of 9 which corresponds to 9 = 0. We will prove that 
dQtop is also a level curve of 9. In other words, computing the value of 9 at the endpoint of 
a level curve emanating from dQt, ase is independent of the level curve chosen. The proof is 
an application of the first Green identity (see Proposition II. 6p . 

Proposition 2.13. The curve dQ top is a level curve of 9 in Q s ut- 

Proof. Let L\ and L 2 be any two level curves of g which start at <9<2base and have their 
endpoints X\ and x 2 on dQ top , respectively. Without loss of generality, assume that the 
g- value of L 2 is bigger than the g- value of L\. We must show that 



First, add vertices of type II according to the procedure defined in Subsection II .3| so that 
the first Green identity, Proposition 11.61 ma Y be applied to the network which is defined on 
the quadrilateral determined by Li,dQ top , L 2 and <9<2base- Let us denote this quadrilateral 

by Q(Li,l 2 )- 

Let h = 1 be the constant function defined in Q(l 1 ,l 2 )- The assertion of Proposition 11.61 
applied with the functions h and g on the induced network, yields 



(2.14) 



0{xi) 



9{x 2 ). 



(2.15) 




However, since g is harmonic in A and in particular along slit (.A), we have that 





Therefore, we must have that 



(2.17) 
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Since g is harmonic in A, it follows easily that 
(2.18) £^ t^(Qe 2 m)( x ) + J ^(Q( Ll ,L 2) )(x) = 0. 

Equation fl2TT4j) now follows from (l2~T7j) . ( 12~18|) and the definition of 9 (Definition UM>- 

□ 

As a consequence of this proposition we can now make 
Definition 2.19. The period of 9 is defined by 

(2.20) period^) = [ ^(Q* 2 ,zko)(«), where m E {0, . . . , k}. 

Ju£L(v m ) an 

We continue the study of the level curves of 9. Note that the maximum principle, and its 
definition, 9 is monotone increasing along level curves of g. This property will now be used 
in the following 

Proposition 2.21. Each level curve of 9 has no endpoint in the interior of Q s ut, is simple, 
and joins Ei to E 2 (9Qbottom t° dQ top in the quadrilateral case). Furthermore, any two level 
curves of 9 are disjoint. 

Proof. Suppose that a level curve of 9 which starts at s E E 2 has an endpoint £ in T E 1~^ 2 \ 
where T lies in the interior of A. Let [s,£] be the intersection of this level curve with the 
interior of A. Let denote the level curve of g that passes through £. Since the level curves 
of g foliate A, there exists a level curve of g, which is as close as we wish to Lg, and 
such that its intersection with [s, £] is empty. Since 9 is monotone increasing and continuous 
along L^, it assmes all values between and period^). Hence, it will assume the value #(£). 
This shows that no level curve of 9 can have an interior endpoint. 

Assume that one of the level curves of 9 is not simple. Let T> be any domain which is 
bounded by it. Since the level curves of g foliate the annulus, one of these intersects the 
boundary of T> in at least two points. The monotonicity of 9 along the level curves of g 
renders this impossible. 

Assume that there exists a level curve of 9, L(9) which does not join E\ to E 2 . By 
construction, each level curve of 9 does not have an endpoint inside A and its intersection 
with each 2-cell is a segment (or a point). Hence, both endpoints of L{9) must lie on Ei or 
on E 2 . Without loss of generality, assume that both endpoints are on E\. Hence, there must 
be a level curve of g that intersects L(9) in at least two points. This is impossible again. 

The fact that level curves of 9 that correspond to the same value may not intersect each 
other follows from similar arguments to those appearing the first part of the proof. 

□ 

Of special importance is the interaction between the level curves of 9 and the level curves 
of g. The following proposition will show that, from a topological point of view, the union 
of the two families of level curves resemble a planar coordinate system. This proposition 
is the topological prerequisite for the proof of Theorem I0.3[ which will appear in the next 
subsection. 
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Proposition 2.22. The intersection number between any level curve of 9 and any level curve 
of g is equal to 1. 

Proof. It easily follows for the proof of Proposition 12.211 that the intersection number of any 
level curve of g with any level curve of 9 is at most equal to 1. Since both families of level 
curves foliate <2 s iit, the interesection number is equal to 1. 

□ 

2.2. Proof of Theorem I0.3L We now turn to the proof of Theorem 10.31 As required, we 
will construct 1Z and then prove the existence of conductance constants defined on llS 1 ' so 
that both g and 9 are harmonic. It will follow from the construction that Equation (10. 2 j) (see 
Definition 10. ip holds. 

Proof. Each vertex in or a vertex of type I in Q s \ it belongs to one and only one of the 
level curves of 9. Let 

(2.23) V = {L(9 ),... J L(9 m )} 

be the set of level curves of 9 which contain all of the vertices mentioned above. Note that 
by definition L(9 ) = <9Qbase is the level set of 9 which corresponds to the value 0, and that 
by Proposition I2.13[ L(9 m ) is the level curve of 9 which corresponds to the value period(#). 

Let 1Z be the cellular decomposition induced by the union of V and C It follows from 
Proposition I2.22I that each 2-cell in 1Z is a quadrilateral. Note that IZ^ contains more 
vertices than the ones used to define the families C, V. To end the proof of the theorem we 
now prove 

Proposition 2.24. There exists a choice of conductance constants on so that both g 
and 9 are harmonic. 

Proof. The proof is a direct consequence of Proposition 12.221 It implies that each interior 
vertex in TZ^°\ has exactly two edges that are parts of a level curve of g and exactly two 
edges that are parts of a level curve of 9, emanating from it. 

Let v(i,j) denote the vertex which is the intersection of L(vi) and L(9j), with i = 
1, . . . , (k — 2), j — 0, . . . , (m — 2). We define, starting with i = 1 and j = 1, 

(2.25) c(Uij, Vy+i) and C^y^y-i) 

to be any two positive constants such that 

(2.26) c(v itj , v iJ+1 ) (9(vij) - 9(v iJ+1 )) + c{v id , (#0m) ~ = 0. 

Choosing such constants is possible since 9 is monotone increasing in the counter clockwise 
direction along L(vi). 
Similarly, we define 

(2.27) civi-x^ViJ) and c(vij,v i+1J ) 
to be any two positive constants such that 

(2.28) c(vi-ij,Vij) (g(vij) - g(v^ ld )) + c(v i+1:j , v i:j ) (g{v^) - g(v i+1J )) = 0. 

Choosing such constants is possible because g is monotone increasing along every level set 
of 9. 
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Each level set of g or of 9 is homeomorphic to a closed interval. Hence, the choices made 
above can successively be continued along the level curves of g and 9 until the second to last 
edges, [t>i, m _i, v i>m ],i = 1, ...,(k- 2) and [u fc -i, m -i, Vfc_i, m ], m = 1, . . . , (m - 2), respectively. 
For the very last edges, compatible conductance constants are determined, by the all the 
choices that were made before. 

A simple computation now shows that with respect to these choices of conductance con- 
stants on TZ {1 \ both g and 9 are harmonic. 



Proposition 12.241 



Finally, since each edge e E is part of either a level set of g or a level set of 9, it is 
evident that Equation (10. 2p holds. This ends the proof of the theorem. 

□ 

Remark 2.29. It follows from Proposition ^. 13l that the net R can be extended to a rectangular 
net on A. 

It follows from the proof that there is a fair amount of flexibility in the choice conductance 
constants. For the application of this paper, a particular choice hereby will be made-one 
that keeps the L\ norm of the normal derivative of g with respect to unchanged. 

Lemma 2.30. There is a choice of conductance constants on such that both g and 9 
are harmonic, and furthermore 

(2.31) / d J-M = I 



T io) nEl dn Jn^nEi dn 

Proof. In light of the previous proposition, it suffices to show that one can choose conductance 
constants 

(2.32) c(vk-ij, Vk,j), j = 1, • • • , m - 1, 

so that Equation (12.311) holds, and then proceed as before. 
Suppose that 

(2-33) = r. 

It then easily follows by setting (for example) 



J'. 



dg 

TWnEi dn 



V 



(2.34) c{v k -u,v ktj ) 

[m — 1) [k — r) 

that Equation (I2.3ip holds. 



m 



□ 



Remark 2.35. The assertion of this proposition means that computing period(#) (see Def- 
inition I2.19p . with respect to endowed with the conductance constants given in the 
Dirichlet data for g or with respect to IZ^ ' endowed with the conductance constants chosen 
above, yields the same number. 
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2.2.1. Viewing 9 from a PDE perspective. The term "conjugate" associated with 9 is moti- 
vated by the following properties of 9 with respect to 1Z and Q s iit : 

• 9 is harmonic, 

• ^ISbasc = °> 9 1 Stop = period(#), and 

de 

Hence, 9 satisfies the combinatorial analogues of the analytical properties of the polar 
angle function v (r, 0) = in the complex plane, which is known to be, when it is single- 
value defined, the harmonic conjugate function of v (r, 0) = r. 

2.2.2. Related work. Our definition of the harmonic conjugate function is motivated by the 
fact that, in the smooth category, a conformal map is determined by its real and imaginary 
parts, which are known to be harmonic conjugates. The search for discrete approximation of 
conformal maps has a long and rich history. We refer to [21] and [T2J Section 2] as excellent 
recent accounts. 

Convincing evidence that a discrete scheme may converge to the Riemann mapping is 
supplied, for example, by the solution to a conjecture suggested by Thurston and based 
on the Circle Packing Theorem by Koebe-Andreev-Thurston. Thurston conjectured that a 
certain discrete scheme will converge to the Riemann mapping. The conjecture, which was 
proved in 1987 by Rodin and Sullivan (see [25]), provides a refreshing and geometric view 
on Riemann's Mapping Theorem. 

We should also mention that a search for a combinatorial Hodge star operator has recently 
gained much attention and is closely related to the construction of a harmonic conjugate 
function. We refer the reader to [21] and to [25] for further details and examples for such 
operators.. 

2.3. The pair-flux length. In this subsection, we will define a notion of length for the level 
curves of g and of 9. To give some perspective, recall that in [T5], Duffin defined a metric to 
be a function r : E — > [0, oo]. More recently, in [TU], Cannon defined a discrete metric to be 
a function p : V — > [0, oo). The length of a path is then given by integrating r, p along it, 
respectively. 

In [18j Definition 1.9], we defined a metric (in Cannon's sense) which utilized g, the solution 
of the boundary value problem, alone. 

Motivated by the planar Riemannian case (see Equation (I2.4ip ). we will define new notions 
of metric and length for level curves of g in Q s i it and Q, respectively. These notions will 
incorporate both g and its conjugate function, 9. 

Definition 2.36. With the notation of the previous sections we define the following. 

(1) For e = [e~,e + ], let -0(e) = e~ be the map which associates to an edge its initial 
vertex. The pair-flux weight of e is defined by 

(2.37) 

2n 2n 2ir 2n 

P(6) = ^^^^m 9 ^ 1 ^ 1 = p^rio^^p^rio^^ 
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(2) Let L be any path in 1Z\ then its length with respect to the pair-flux weight is given 
by integrating p along it, 



(2.38) Length(L) = / p(e). 

In the applications of this paper, we will use the pair-flux weight to provide a notion of 
length to level curves of g. Thus, by the assertion of Proposition 12. 13[ we may now deduce 

Corollary 2.39. Let L(v) be a closed level curve of g, and let < m = g\i(v) < k; then we 
have 

2ir 

(2.40) Length(L(y)) = 2vrexp ( — m). 

v penod(6 l ) ' 

The definition of the pair-flux length is one of the new advances of this paper. Whereas 
in [18~| [T9| [20] other notions of lengths utilizing only the solution g were introduced, the 
pair-flux length incorporates the pair {g,0}. This appealing feature is also the case in the 
smooth category, i.e., for z = rexp (i9) in the complex plane, we have 

(2.41) dz = r exp {iff)d0 + exp (i6)dr. 

We will now provide a notion of length to the level curves of 9 in Q B m,- Keeping the 
analogy with the planar Riemannian case, the restriction of the Euclidean length element to 
level curves of the function v(r, <fi) = <fro, has the form 

(2.42) \dz\ = dr. 

Definition 2.43. Let L{9) = (v , . . . ,Vf.) be a level curve of 9 with t> G E 2 and Vk G E±, 
then its length is given by 

(2.44) Length(L(^)) = exp(g(v k )) - exp(^(t; )) = exp(A;) - 1. 

Remark 2.45. It is a consequence of part (2) of Proposition 12.221 that any two level curves 
of 9 have the same length. 

3. The cases of an annulus and an annulus with one singular boundary 

component 

3.1. The case of an annulus. The setting of this subsection is a special case of the one 
described in Definition II .41 Let 1Z be the cellular decomposition into quadrilaterals associated 
with the orthogonal pair {g, 9} which was constructed in the proof of Theorem 10. 3[ 

We use the term measure on the space of quadrilaterals in IZ^ to denote a non-negative 
set function defined on IZ^ . An example of such, which will be used in Theorem I0.4[ is 
provided in 

Definition 3.1. For any R G TZ^\ let Rtop, Rbasc be the top and base boundaries of R, 
respectively. Let t G -R^p and b G R^ sc be any two vertices. Then we let 

/o r\\ / D \ 1/ 2/ 27T / 1 \ \ 2/ 2n , y ^2nd9(R hase ) 

(3 ' 2) V{R) = 2^ { ^d(9) 9{t)) " eXP ( p^rio^ (&)) ) period(fl) ' 
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Remark 3.3. By construction, all the vertices in R top (i?base) have the same g values and 
d9(Rbase) — d9{R top ). 

We now turn to the 

Proof of Theorem I0.4L Recall (see the discussion preceding Proposition 12.241) that the 
vertices in IZ^ are comprised of all the intersections of the level curves of g (the family C) 
and the level curves of 9 (the family V). The vertex will denote the unique vertex 
determined by the intersection of L(v{) and L(9j) (the existence and uniqueness of this 
intersection are consequences of Proposition 12. 22) ) . 

The conjugate function 9 is single- valued on Q s u t and multi- valued with a period which is 
equal to period(#), when extended to A. This means that 

(3.4) 9( Zl ) = 9(z )+penod(9), 

whenever z\ G L(zq) is obtained from zq G slit (A) by traveling one full cycle along the level 
curve L(zq). Hence, the discrete function 

(3.5) period(g) (9(v)+i9(v)), v G A n fc<°> 
has period 2iri when defined on A. Therefore, 

(3.6) exp(— (g(v)+i9(v))), v G A D 

is single- valued on A. 

We now turn to the construction of the tiling T. 

The tiling T of Sj± is determined by all the intersections of the family of concentric circles, 
C, defined by 

(3.7) n = exp ( periQd( g ) 9(vi)) , for i = 0, . . . , k, 
with the family of radial lines T, defined by 

(3.8) t^-JlL-Oj, for J = 0,...,m, 

where each annular shell in the tiling is uniquely defined by four vertices that lie on two 
consecutive members of the families above. 

Let f-ji be the homeomorphism which sends the quadrilateral TZ determined by the coun- 
terclockwise oriented ordered set of vertices 

(3.9) (i + l,j), (i + 1, j + 1), + 1)} 

for % — 0, . . . , k and j = 0, . . . , m, onto the annular shell Tr determined by the counterclock- 
wise oriented ordered set of vertices 

(3.10) {rjexp (z^-), r i+1 exp (z^), r i+1 exp (i<fr j+ i), r { exp (ifa+i)} 
and that preserves the order of the vertices. 

By the definitions of v, fx and Tr, we have for all R G that 

(3.11) v{R) = fx{T R ). 
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Let 

(3.12) /= |J f R . 

Then it is clear from (13. 9 p and ( 13 . 1 Of) that / maps any edge in TZS 1 ' D dVL homeomorphically 
onto an arc in dS^. Therefore, / is indeed boundary preserving. 




Figure 3.13. Constructing one annular shell. 



It remains to prove that the map / = Ur/h assembled from the individual maps defined 
above, which is clearly a homeomorphism onto its image, is onto S^. To this end, observe 
that by the maximum principle, the map / is into S 1 ^. 

It is clear from the construction of 1Z that any two quadrilaterals in TZ^ have disjoint 
interiors and that their intersection is either a single vertex or a common edge. Also recall 
that by definition each quadrilateral has its top and bottom edges situated on two successive 
level curves in C, and its right and left edges situated on two successive level curves in V. 
Since the union of the quadrilaterals in tile A, the total //-measure of their union, which 
we define to be v(A), satisfies 

(3.14) u{A) = u{ (J R)= 

Starting from the quadrilaterals that lie between L(vq) and L(v\), we sum the ^-measure 
of all the quadrilaterals in the layer defined in between successive level curves of g, until we 
reach L{v^). An easy computation employing Definition I3. II now shows that 



(3.16) 
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By the construction of the annular shells and the definition of the map /, each quadrilateral 
R is mapped onto a unique annular shell T R , no two different quadrilaterals are mapped onto 
the same annular shell, and the collection of their images tiles a subset of S^. 

Hence, by applying the above paragraph, (13.141) . (13. lip . (I3.15P and the definition of S^, 
we obtain that 

(3.16) u(A)= J2 KTr)=I*( (J T r ) = im(S a ). 

Hence, there are no gaps nor overlaps in the tiling of 1S4. This concludes the proof of the 
theorem. 

□ 

Remark 3.17. The proof further shows that in fact level curves of g are mapped homemorphi- 
cally onto the level curves u(r, <p) = r*, and level curves of 9 are mapped homeomorphically 
onto the level curves v(r, 0) = <f>j. 

3.1.1. Relation to works by Schramm and Cannon-Floyd-Parry. It seems imperative to relate 
this theorem to Theorem 1.3 in [27] and Theorem 3.0.1 in [IT] . While Schramm, and by 
Cannon, Floyd, and Parry used discrete extremal lengths arguments in their proofs, the 
arguments as well as the results are a bit different. Schramm's proof seems to work for a 
quadrilateral but not directly for an annulus. The methods of Cannon, Floyd and Parry work 
for both a quadrilateral and an annulus. Furthermore, Schramm's input is a triangulation 
with a contact graph that will (more or less) be preserved. The input for Cannon, Floyd and 
Parry is more flexible. They consider a combinatorial of a topological quadrilateral (annulus) 
by topological disks. We refer the reader to the papers above for details. Upon applying 
a Dirichelt-Neumann boundary value problem, our proof will work for the quadrilateral as 
well. 

It is showed in [27J page 117] that if one attempts to use the combinatorics of the hexag- 
onal lattice, square tilings (as provided by Schramm's method) cannot be used as discrete 
approximations for the Riemann mapping. There is still much effort by Cannon, Floyd and 
Parry to provide sufficient conditions under which their method will converge to a conformal 
map. 

While our proof of Theorem 10.41 does not use the machinery of extremal length arguments, 
it is worth recalling that in the smooth category there are celebrated connections between 
boundary value problems and extremal length (see for example [U Theorem 4.5]). 

The common theme of our methods and those of Cannon, Floyd and Parry in [TTJ is 
the construction of a new coordinate system on a topological annulus. As stated in the 
introduction, this idea goes back to Riemann. 

3.2. The case of an annulus with one singular boundary component. In this subsec- 
tion we will generalize Theorem 10.41 by providing a geometric model for an annulus with one 
singular boundary component. The singular boundary component is of a special type. It is 
determined by the topological structure of a singular level curve of the solution of a Dirichlet 
boundary value problem imposed on a planar embedded m-connected domain, where m > 1. 
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We start with two definitions; the first one appeared in [THJ page 9]. 

Definition 3.18. A generalized bouquet circles will denote a union of bouquets of piecewise- 
linear circles where the intersection of any two circles is at most a vertex. Moreover, all such 
tangencies are required to be exterior, i.e., no circle is contained in the interior of the bounded 
component of another. 

Recall that Theorem 2.15, which was proved in [TS], asserts the following. 

Theorem 3.19 (The topology of a level curve). Let L be a connected level curve for g. Then 
each connected component of L is a generalized bouquet of circles. 

It is convenient to present the singular boundary component as a quotient space. In the 
following definition, a circle will mean either a round circle or a piecewise linear circle. 
To this end we make 

Definition 3.20. An embedded planar circle with finitely many distinguished points on it, 
will be called a labeled circle. If in addition, equivalence relations among these points is 
given, so that the quotient of the labeled circle is a generalized bouquet of circles, then we 
call the quotient a labeled bouquet, and the labeled circle will be called good. 

Remark 3.21. Note that if a labeled round bouquet, i.e., one which consists of only round 
circles, has more than two round circles tangent at one point, it will no longer embed in R 2 . 

We will now define the object of study in this subsection. By a generalized singular annulus, 
^4 s j n g, we will mean a subset of the plane, whose interior is homeomorphic to the interior of 
an annulus, and whose boundary has two components: one of which is homeomorphic to S 1 
and the other is a generalized bouquet of circles. The subscript denotes the set of tangency 
points in the generalized bouquet of circles. Let us also assume that a cellular decomposition 
T of (-4 s i n g, ^-sing) ^ s provided, where each 2-cell is either a triangle or a quadrilateral. 

Topologically, A s [ U g may be presented as the quotient of a planar annulus A, where 
dA = Si U £?2) and E2 is a good labeled circle (see Definition I3.20p . Henceforth, we will 
let it denote the quotient map. We will let E2 denote the singular boundary component of 
8 A ■ 

Note that the cellular decomposition T can be lifted to a cellular decomposition T of 
(A,8A), where each 1-cell, 2-cell in T, respectively, is the unique pre-image, under 7T" 1 , of 
a unique 1-cell, 2-cell in T, respectively. The difference between the two cellular decompo- 
sitions manifests in the addition (in comparision to E2) of vertices in E2. Specifically, for 
each vertex v in the singular part of 8A s - m g, there are m(v) vertices in E2, where m(v) is 
the number of circles that are tangent at v. 

We will now apply Theorem 10.31 and Theorem 10.41 to (A, 8 A, T). In the following proposi- 
tion, recall that the existence of 1Z is provided by Theorem I0.3[ and that 9 is the conjugate 
harmonic function to g, the solution of the imposed Dirichlet boundary value problem on 
(A,8A,f). 

With the above notation and setting in place, and with £ denoting the set of level curves 
of g as in equation (12.31) . we may now state the main proposition of this subsection. 
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Proposition 3.22 (An annulus with a singular boundary). Let (-4. s i n g, cM- S mg = E\ U E 2 ) 
be a generalized singular annulus endowed with a cellular decomposition T ■ Let k be a 
positive constant and let g be the solution of the Dirichlet boundary value problem defined on 
(A, dA, f). 

Let S A be the concentric Euclidean annulus with its inner and outer radii satisfying 

2ir 

(3.23) {n,r 2 } = {1,2tt Length^))} = {l,2vrexp ( ——- ft)}. 

Then there exist 

(1) a tiling T of S A by annular shells, 

(2) a set denoted by sing^^) consisting of finitely many points which is contained in the 
inner boundary of S A , 

(3) a homeomorphism 

f : (A, dA \ TT-^singGA)), U) (S A , dS A \ sing(SU), T) 

such that f maps the interior of each quadrilateral in IZ^ onto the interior of a 
single annular shell in S A , f preserves the measure of each quadrilateral, i.e., 

u(R) = (i(f(R)),ioi all R G 7Z {2) , 

and f is boundary preserving. 

Proof. The proof is a straightforward modification of the non-singular boundary case. Let 
1Z be the rectangular net constructed in Theorem IU .31 Let 9 be the conjugate harmonic 
function constructed on (A, dA, T), and let / be the map constructed in Theorem 10.41 Let 
T be the tiling of S A provided by Theorem 10.41 

For each ti G sing(^4.) C E 2 , % = 1, . . .p, there are precisely m(tj) vertices on E 2 in the 
equivalence class corresponding to ti. Let 

(3-24) V{U) = {L{6) Kl ,. . . , L(6) t . m(W }, i = l,...p 

be the level curves of 9 that have one of their endpoints at these vertices. With this notation, 
and since the level curves of 9 are parallel, it follows that 

p 

(3-25) V sing( ^ = |J V ^) 

i=l 

comprises of all the level curves of 9 that have an endpoint in the pre- image of sing(^4). 
Set 

(3.26) sing(^) = f(E 2 )f]f((jV(t t )), 

i=i 

then sing(S'^4) is the image under / of all the vertices in the pre-image of sing(^4). Further- 
more, recall that f0^sing(s A )^ ^ s a se ^ °^ racua l arcs ^ n which are parts of TZ. 

To finish proving the statement in (3), note any quadrilateral in 7Z^ whose vertices are 
disjoint from sing(S , ^i) is mapped homemorphically onto a shell in S A . Since by construction 
the image of 7r _1 (sing(^4)) is precisely sing^^), it follows that / will map the interior of 
each one of the rest of the quadrilaterals homemorphically onto the appropriate shell, with 
punctures at the corresponding vertices. This ends the proof of the proposition. □ 
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A geometric model to («4 s j n g, dA s [ n g, T) is now easy to provide since the first part of 
(3) in the proposition above allows us to label the vertices in sing^^) isomorphically to the 
labeling of the vertices in sing(^4). We will keep denoting by 7r the quotient map which is 
thereafter induced on Sjs,. Such a quotient annulus will be called a generalized Euclidean 
annulus and will be denoted by Cj±. The proof of the following corollary is straightforward. 

Corollary 3.27. With the assumptions of Proposition \3. 22\ and with = S^/n, there 
exist 

(1) a tiling T of by annular shells, and 

(2) a homeomorphism 

f '■ (-^sing' ^sing'^) ~~ > \^AjdC^,T), 

such that /(sing(^4)) = sing(S'^)/7r ; / maps the interior of each quadrilateral in 
onto the interior of a single annular shell in Cj±, f preserves the measure of each 
quadrilateral, i.e., 

u(R) = n(f(R)), for all R G TZ {2 \ 
and f is boundary preserving. 




Figure 3.28. An example of the map ir. 



In the next section, we will work with a general m-connected planar domain that will 
be cut along singular level curves of a Dirichlet boundary value problem imposed on it. 
Euclidean cylinders and Euclidean cylinders with one singular boundary component will be 
utilized in order to allow gluing of two components of a singular level curve. To this end, 
recall that a conformal homeomorphism, from a concentric annulus to a Euclidean cylinder 
of radius equals to 1 and height equals to log (a/b), is defined by 

(3.29) F(pcos(#),psin(#)) = (cos (0), sin (0), log (p)), a < p < b, < 9 < 2vr, 
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where (p, 6) denote polar coordinates in the plane. 

It easily follows that the image of an annular shell under the map F is a Euclidean 
rectangle. We will abuse notation and will keep the same notation for C4 and its image 
under the mapping F. 

We now define a variation of the measure v (see equation (j3.2p ) in order to adjust our 
statements to working with such cylinders. 

Definition 3.30. For any R £ 7o 2 \ let Rt op ,Rbase be the top and base boundaries of R, 
respectively. Let t £ R^ p and b £ R^ sc be any two vertices. Then let 

(3 31) \(R) - 2rrde ( Rhas ^ fog 9(t) 

period(6 l ) g{b) 

By applying the map F and the measure A, we may state Theorem 10 A\ Proposition 13.221 
and Corollary 13.271 in the language of Euclidean cylinders. We end this subsection by sum- 
marizing that in the following remark which will be applied in the next section. 

Remark 3.32. Under the assumptions of Theorem 10.41 Proposition 13.221 and Corollary 13. 27[ 
all the assertions hold if one replaces S^, generalized by a Euclidean cylinder, generalized 
Euclidean cylinder, respectively; / by ir o F o f and an annular shell by its image under F, 
F o 7r, respectively, and the measure v by the measure A. 

4. PLANAR DOMAINS OF HIGHER CONNECTIVITY 

In this section we will generalize Theorem 10.41 to the case of bounded planar domains of 
higher connectivity. We will start by recalling an important property of the level curves of 
the solution of the Dirichlet boundary value problem (see Definition II ,4p . This property will 
be essential in the proof of Theorem l4.41 In the course of the proof, we will need to know that 
there is a singular level curve which encloses all of the interior components of <9fi, where Q 
is the given domain. This level curve is the one along which we will cut the domain. We will 
keep splitting along a sequence of singular level curves in subdomains of smaller connectivity 
until the remaining pieces are annuli or generalized singular annuli. Once this is achieved, 
we will provide a gluing scheme in order to fit the pieces together in a geometric way. 

Before stating the second main theorem of this paper, we need to recall a definition and 
a proposition. Following [2] and [231 Section 3], consider the number of sign changes in the 
sequence {/(wi) - f(v), f(w 2 ) - f(y), . . . , f(w k ) - f(v), f(w 1 ) - f(v)}, which is denoted by 
SgCj(f). The index of v is then defined by 



(4.1) lnd f (v) 



Sgc f {v) 



Definition 4.2. A vertex whose index is different from zero will be called singular; otherwise 
the vertex is regular. A level set which contains at least one singular vertex will be called 
singular; otherwise the level set will be called regular. 

The following proposition appeared as Proposition 2.28 in [18]. 

Proposition 4.3. There exists a unique singular level curve which contains, in the interior 
of the domain it bounds, all of the inner boundary components of dQ. 
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Such a curve will be called the maximal singular level curve with respect to f2. Recall that 
the notion of an interior domain was discussed in Subsection 11.31 

Throughout this paper, we will not distinguish between a Euclidean rectangle and its 
image under an isometry. Recall (see the end of Subsection I0.3j) that a singular flat, genus 
zero compact surface with m > 2 boundary components with conical singularities is called 
a ladder of singular pairs of pants. 

We now state and prove a generalization of Theorem 10.41 

Theorem 4.4 (A Dirichlet model for an m-connected domain). Let (Q, dfl = E\ U E 2 , T) be 

a bounded, m-connected, planar domain with E 2 = E 2 UE 2 . . . U E™' 1 . Let g be the solution 
of the Dirichlet boundary value problem defined on (Q,dQ,T). Then there exist 

(1) a finite decomposition A = UjA; of Q, where for all i, Ai is either an annulus or an 
annulus with one singular boundary component; 

(2) for all i, a cellular decomposition TZA^f Ai where each 1-cell is a simple quadrilateral; 

(3) for all i, a finite measure A, defined on TZaJ and 

(4) a ladder of singular pairs of pants Sn with m boundary components, such that 

(a) the lengths of the m boundary components of Sq are determined by the Dirichlet 
data, 

(b) there exists a decomposition of Sq = UjC^, where each Gj^ is either a Eu- 
clidean cylinder or a generalized Euclidean cylinder, equipped with a tiling T; t by 
Euclidean rectangles endowed with Lebesgue measure; and 

(c) a homeomorphism 

/:(n,9fi,u i ^ i )H>(5 n ,9n,u i ^ Ai ), 

such that f maps each Ai homeomorphically onto a corresponding C^, and 
each quadrilateral in TZa% onto a rectangle in while preserving its measure. 
Furthermore, f is boundary preserving (as explained in Theorem \0.J$ . 

Proof. The first part of the proof is based on a splitting scheme along a family of singular 
level curves of g which will be proven to terminate after finitely many steps. We will describe 
in detail the first two steps of the scheme, explain why it terminates, and leave the "indices" 
bookkeeping required in the formal inductive step to the reader. The outcome of the first part 
of the proof is a scheme describing a splitting of the top domain, Q, to simpler components, 
annuli and singular annuli. 

The complement of L(Q), the maximal singular curve in Q, has m-connected components, 
all of which, due to Proposition 14.31 have connectivity which is smaller than m — 1, or are 
annuli, or generalized singular annuli. By the maximum principle, one of these components 
has all of its vertices with g- values that are greater than the g- value along L(Q). In Subsec- 
tion [L3l such a domain was denoted by 2 (L(Q)) and was called an exterior domain. Its 
boundary consists of E\ and L(Q). It follows from Proposition 14.31 and Theorem 13.191 that 
it is a generalized singular annulus which will be denoted by A(Ei, L(Q)). 

Let the full list of components be enumerated as 

(4.5) Cd = {CC 1A (L(ty), CC 1)2 (L(n)), . . . , CC ljm (L(ty) = A(E U L(Q))}. 
By definition, for each j — 1, . . . ,m — 1, the g- value on the boundary component 

(4.6) d ld = dCCu(L(n)) n L(Q) 
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is the constant which equals the g-value on L(Q). The other components of dCCij(L(Q)), 
j = l,...,m — 1, are kept at g-values equal to 0. Hence, we now impose a Dirichlet 
boundary value problem with these values on each element in the list CC\ \ A(E\, L(Q)). 
On A(Ei, L(Q)), the induced Dirichlet boundary value problem is determined by the value of 
g restricted to Ei and the value of g restricted to L(Q). Note that imposing these boundary 
value problems in general will require introducing vertices of type I and of type II and 
changing conductance constants along new edges, as described in Subsection 11.31 These 
modifications are done in such a way that the restriction of the original g solves the new 
boundary value problems. 

For j = 1, ...,m — 1, let kj denote the connectivity of CCij(L(Q)). We now repeat 
the procedure described in the first paragraph of the proof in each one of the connected 
components CCij(L{Q)), kj — 2 times, for j = 1, ... ,m — 1, excluding those indices that 
correspond to annuli. 

We will now describe the second step of the splitting scheme. For each j £ {1, . . . , m — 1} 
whose corresponding component is not an annulus, a maximal singular level curve Lj(CC\j) = 
L(CCij(L(Q))) with respect to the component CCxj(L(JY)) is chosen. This is possible be- 
cause at the end of the previous step we imposed a Dirichlet boundary value problem on 
each one of these domains. Hence, the assertion of Proposition 14.31 and Theorem 13.191 may 
be applied to these domains as well. 

Therefore, a new list consisting of connected components of the complement of Lj(CC\J) 
in CC\j(L(Q)), of cardinality at most m — 1, 

(4.7) CCij = {CCi,j,i{Lj{CCij))i CC\j2,{Lj{CCiJ)), . . . , CCi^m-^L^CCij))}, 

j as chosen above is generated. We will let the last element in this list denote the exterior 
domain to Lj(CC\j) in CCij. It is, as in the first step of the scheme, a generalized singular 
annulus denoted by A(dij, Lj(CCij)). The other components have connectivity which is 
smaller than (m — 2), or are annuli. Note that in this step the exterior domain from the first 
step in the scheme, CCi im (L(Q)) = A(Ei, L(Q)), is left without any further splitting, since 
it is a generalized singular annulus. 

By definition, for each j chosen as above, and each % = 1, ... m — 2, the g-value on the 
boundary component 

(4.9) d ld>i = i)C(\. l JL J (CC l . J )) n LjiCCxj) 

is the constant which equals the g- value on Lj(CCij). The other components of dCCij t i(L(Q)) 
are kept at ^-values equal to 0. Hence, we now impose a Dirichlet boundary value problem 
with the above values on each connected component in the list CCij \ A(dij, Lj(CCij)). 
On A(d\j, Lj(CCij)) the induced Dirichlet boundary value problem is determined by the 
value of g restricted to d\j and the value of g restricted to Lj(CCij). Addition of vertices 
of type I and II, and modifications of conductance constants will be applied in this step as 
in the previous one. 

It follows that in each step of the splitting scheme either a domain with a fewer boundary 
components than the one that was split, or a singular annulus, or an annulus, is obtained. 
Hence, the connectivity level of each connected component after the split is either the con- 
stant number two, the constant number three, or it decreases. Therefore, the splitting scheme 
will terminate once all the obtained components have connectivity which equals to two or 
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Figure 4.8. An example of a splitting scheme. 



three, i.e., when the union of all the final generated lists is a list of lists, each containing 
only annuli and generalized singular annuli. 

We now turn to the second part of the proof. Here, we will show that it is possible to 
reverse the splitting scheme, i.e., starting at the final lists generated in the splitting scheme 
up to the first one, CC\, we will glue the pieces in a geometric way; that is, so that the 
lengths of glued boundary components are equal. It is in this step where the pair-flux length 
will be used. 

By the structure of the lists obtained in the first part of the proof, it is sufficient to show 
how to 

(1) glue in a geometric way elements in a list, say CCi^...^, that contains only annuli 
and a generalized singular annuli so as to form a ladder of singular pair of pants 
denoted by <Sij,. an d 

(2) glue in geometric way «Si lJ ,... ; fe to the singular boundary component of the generalized 
singular annulus in the list from which CCi^...^ was formed. 

Note that, if all the elements in the list CCi^...^ are annuli, we can apply case (2), 
since we will map each annulus via Theorem 10.41 to a Euclidean cylinder (by first applying 
Theorem 10.41 and then Remark 13.321) . In order to ease the notation, let us show steps (1) 
and (2) for the lists produced in the first part of the proof. 

Fix a j G {1, . . . , m — 1}, and for i = 1, . . . , m — 2, we apply Theorem 10.41 to A% = 
CCijj, which by the assumption of step (1) is an annulus. This yields a collection of 
concentric Euclidean annuli {Si, . . . , S m -2} where, in the notation of Theorem 10 A\ Si = S^. 
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Furthermore, for each i = 1, . . . ,m — 2, we have 

(4.10) {r\y 2 } = {MTrLength^)^^))} = {l,27rexp (-^—. g (d lJti ))}, 

where 6i is the period of the conjugate function to g, the solution of the Dirichlet boundary 
value problem defined on CCi^j. 

We now apply the map F defined in equation ( I3.29P to obtain a corresponding sequence 
of Euclidean cylinders {C\, . . . , C m _2}- All of these have radii equal to the constant 1, and 
their heights are given respectively by 

(4.11) hi = log (27rLength (gA) (L(a 1 , J -, i ))) = log27r + 

Recall that the last component in the list CC\j is the generalized singular annulus 
A(d\j, Lj(CCij)). Let 9xj be the conjugate harmonic function to g, the solution of the 
Dirichlet boundary value problem induced on it (see the paragraph preceding Proposi- 
tion [3722]) . Then, following Remark 13.321 we now map it to a generalized Euclidean cylinder, 

Measuring lengths of the image of the boundary components of A(dij, Lj(CCij)) (using g 
and 6i j)i which are in one to one correspondence with the singular boundary of this cylinder, 
yields the sequence of radii 
(4.12) 

1 l-cng.1,,,07,.,,) 1 1 ex P ( periS( g<J ) ^^))) f M n 

1 2tt Lengthy. (d lsj ) 2tt period^) exp (—^- )9 (d ld )) J ee9l ^ 

for % = 1, . . . , m — 2, where the expression on the righthand side is based on equation ( 12. 38ft 
and equation (12.401) . The sequence Ri, i — 1, . . . ,m — 2 comprises the sequence of lengths of 
the round circles in the singular boundary component of C^(dij, Lj(CCij)). 

Let 

(4.13) ^ = f Ri = R iZ 

be the conformal homeomorphism acting on the Euclidean cylinder d = it o F o /(CCij^), 
where z is the standard complex parameter on CC\^^. Hence, we may glue ji[C\) along 
one of its boundary components to the corresponding component in the singular boundary 
component of -^(Cdj), so that the length of the two boundaries are the same. 

This establishes step (1). We will now show how to establish step (2). 

The completion of step (1) yields a ladder of singular pairs of pants which we will denote 
by Si j. This ladder has (m — 1) components; one corresponds to d\j and the others com- 
prise one boundary component of each one of the cylinders fi o d, i = 1, ... m — 2, after 
these cylinders are attached. Recall that d\,j is the intersection of the generalized singu- 
lar annulus A{E\,L(Q)) with the generalized singular annulus in the list CC\j, which is 
A(d 1 , j ,L j (CC lJ )). 

Let g be the solution of the induced Dirichlet boundary value problem on A(E\, L(Q)) and 
let 9 be the harmonic conjugate function to g (see the paragraph preceding Proposition 13. 22p . 
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Let 

(4 ' 14) Tj = Length^-)' 

where lj denotes the length of boundary component that corresponds to dij in Sij, and 
CU(Ei,l(c)) is the generalized Euclidean cylinder constructed for A(Ei, L(Q)) (see Remark [332]). 
After applying a conformal expansion of magnitude Tj to Sxj, it then may be glued along 
this boundary component to C^Ei,L(n)) in such a way that the length of the corresponding 
circle in the round bouquet dC^E lt L(n)) has the same length. This establishes step (2). 

By construction, it is clear that the pair-flux length of the boundary component of Sq 
that corresponds to E\ is equal to 

(415) 2 * exp( i^dW fc) - 

where 9 is the harmonic conjugate to g, the solution of the Dirichelt boundary value problem 
induced on A(Ei, L(Q)) (recall from Definition 11.41 that k is the value of g restricted to E\). 
The lengths of the remaining (m — 1) boundary components of Sq which correspond to the 
m — 1 boundary components of E2 are determined by the process described in the previous 
part of the proof. 

The length of a component in Sq which corresponds to E l 2 , i e {1, . . . , m — 1}, measured 
with respect to the pair-flux metric, is obtained by successively multiplying a sequence of 
ratios of lengths. These ratios are uniquely determined as in equation f)4.12p and equa- 
tion f)4.14p . and present the expansion factor needed in order to match the gluing of a 
(generalized) cylinder to the one which induced it in the splitting process. 

Cone angles are formed whenever more than two cylinders meet at a vertex; viewed in Q, 
this will occur whenever more than two circles in a generalized bouquet meet at a vertex. 
The computation of the cone angles is solely determined by g and T . This analysis first 
appeared in Theorem 0.4 in [18]. Specifically, the cone angle <p(v) at a singular vertex v, 
which is the unique tangency point of n + 1 Euclidean cylinders, staisfies 

(4.16) <f)(v) = 2(n+ 1)tt. 

The proof of the theorem is thus complete, with / defined to be the union of the individual 
maps constructed at each stage. 

□ 

Remark 4.17. There is a technical difficulty in our construction if some pair of adjacent 
vertices of has the same (/-values (the first occurrence is in equation (I4.ip ). This issue 
has already been discussed and resolved in [181 Remark 4.14]. 
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